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Abstract 

We obtain an exact many-body scattering eigenstate in an open quantum dot system. The 
scattering state is not in the form of the Bethe eigenstate in the sense that the wave-number set 
of the incoming plane wave is not conserved during the scattering and many-body bound states 
appear. By using the scattering state, we study the average nonequihbrium current through the 
quantum dot under a finite bias voltage. The current-voltage characteristics that we obtained by 
taking the two-body bound state into account is qualitatively similar to several known results. 
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FIG. 1: A two-electron scattering state which contains incoming plane waves only in the left lead. 

Mesoscopic transport of interacting electrons has attracted much interest recently [2I, I^, 
J]. A remarkable feature of the mesoscopic system is the coherence length greater than the 
sample size. In the standard theory, the electron in the sample is described by the quantum 
mechanics and dissipation is considered to occur only in reservoirs connected to the sample. 
A well-known approach to the electric current across the sample under a finite bias voltage 
is the Landauer formula, although the original one is restricted to the non-interacting case. 



10|. To 



The Green's function is also employed to study the transport property 
discuss the effect of interactions in this framework, however, we would have to resort to a 
perturbation technique, which is generally a hard task. 

In this Letter, we present an exact many-body scattering eigenstate in an open quantum 
dot system and apply the eigenstate to analysis of the nonequilibrium current. The system 
we study is an open interacting resonant-level model (IRLM), which consists of two leads 
of non-interacting spinless electrons that interact with an electron on a quantum dot in 
between the two leads. Each lead is connected to a large reservoir. First, we explicitly 
construct two- and three-electron scattering states, which are free-electronic plane waves 
before scattering and, at the quantum dot, are partially scattered to a many-body bound 
state due to the Coulomb interaction. Second, by using the scattering states, we calculate 
the quantum-mechanical expectation value of the current through the quantum dot in the 
second order of the inverse system length. Third, we study the statistical average of the 
nonequilibrium current for a given finite bias voltage under the assumption that electrons 
are completely thermalized in each reservoir before returning to the lead. 

Our study of the nonequilibrium current with scattering states is a genuine extension of 
the Landauer formula. Our scattering states of the open system are suitable for describing 
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incident electrons thermalized to a free-electron state in each reservoir. Some used the 
Bethe ansatz [ll, Q, Q to study the transport properties of quantum dot systems fl^ . 
where the Landauer formula was formally applied to the quasi-particles in a closed system 
in equilibrium. However, the periodic boundary conditions imposed on the Bethe state are 
clearly different from the conditions adopted for the Landauer formula, the conditions that 
the incident electrons are asymptotically free. Recently, there have been a few attempts to 
study the transport properties with a scattering state in the framework of the Lippmann- 
Schwinger (LS) equation [isl, Q. Our scattering state is shown a solution of the LS equation 
associated with the open IRLM. 

A remarkable point of our solution is the appearance of a many-body bound state in the 



scattering eigenstate. Another many-body bound state given by the Bet 



le ansatz method is 



known to be the ground state of the Anderson model in equilibrium 17|. Our bound state, 
on the other hand, is generated as a result of the scattering of an incident free-electronic 
plane- wave state (Fig. [1]). The interaction around the quantum dot is a necessary condition 
of the appearance of the bound state. The nonequilibrium current is indeed affected by the 
interaction through the bound state. 
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he open IRLM out of equilibrium has been studied with various approaches 
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22|. We express the quantum- mechanical expectation value of the current as a 



series of the inverse system length to consider the average current, while the perturbative 
result 201] gives the average current as a series of the interaction parameter. The qualitative 
behaviour of the current-voltage characteristics that we obtain is similar to the results in 



Refs. 
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2Q|, 



21| . We remark that, in our results, the effect of the interaction appears in 



the quantum-mechanical expectation value, which differs from the result in Ref. 18 1. 
The Hamiltonian of the open IRLM is given by 

dx cl{x)-—Caix)+t{cl{0)d+d^CaiO)) 

L/2 



H J2\ I :~z i dx 

+ eddU + Ucl{Q)c^{Q)dU 



(1) 



where cjj(x) and Ca{x) are creation- and annihilation-operators of the electrons in the lead 
a(= 1, 2), d"^ and d are those in the quantum dot, i = t/ \/2 is the transfer integral between 
each lead and the dot, is the gate energy of the dot and f/(> 0) expresses the Coulomb 
repulsion. The dispersion relation in the leads is linearized in the vicinity of the Fermi 
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energy to be E = v-pk, under the assumption that t, and U are small compared with the 



Fermi energy 
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131 ] ■ For simplicity, we have set f p = 1 in Eq. ([T]). We treat the system 
as an open system in the limit L oo. The lead a is connected infinitely far way to a large 
reservoir characterized by the Fermi distribution with the chemical potential fi^- Our goal 
is to calculate the statistical average of the current 

i = iij2i-r^'{cim-d^cM) (2) 

a 

for the system under a finite bias voltage, fii > /i2- 

We consider the general form of eigenstates. After the transformation Ci/2{x) = [ce{x) ± 
Co{x)) /^/2, the Hamiltonian ([1]) is decomposed into the even and odd parts. Due to the 
relations [H, + N^] = [H, No] = for the number operators N^/o = fdx cl^^{x)ce/o{x) and 
Nd = Sd, the set {Nf. + N^, No} gives a good quantum number. The A^-electron state \N, n) 
in the sector with No = n is expressed in the form 



|iv,n) yrfx^--ci2/V"n^;2/)4(^i)---4(^iv-n)4(i/i)---ct(z/„) 

+ /"rfx^-"-irfy"e(")(x; y)cl{x,) ■ ■ ■ cl{xM-n-i)d^ci{y^) ■ ■ -cKz/n)) |0), (3) 



where we put e^^\x;y) = 0. The functions y) and e^"'\x;y) are antisymmetric with 

respect to the variables {xi} and with respect to {yi}- The one-electron eigenstate |l,r;,; k), 
{n = 0, 1) with the energy eigenvalue E = k is obtained by inserting the eigenfunctions 
g^^\x) = gk{x), e^°^ = Ck or g^^\y) = hk{y) into the general form ([3]), where 

9k{x) = ^e''%e{-x) + ^e{x)), 
V27T e,^ 

hk{x) = ^e'^", efc = ^—tt^, (4) 

V27r V27r k- ed + it^ 



with the step function 0{x). The linear combination \k) = (|1,0; k) + |1, 1; k))/\/2 gives a 
scattering state containing an incoming electron only in the lead 1. If we imposed periodic 
boundary conditions to the leads, the wave number k allowed for the eigenfunction g^^\x) 
would be different from that for g^'^\x). Thus, even in the non-interacting case, the scattering 
state I A;) is inconsistent with the periodic boundary conditions. 

For N = 2, the eigenvalue problem H\2,n) = E\2,n) is cast into a set of the Schrodinger 
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equations: 



}-(^, + ^,)-E)g^'\x^,x,)-^-{6{x^)e^'^\x2)-e^'\x^)6{x2))=0, 
^ +f/5(x) + erf-E)e(°)(x) + 2t(/(°)(x,0) =0, 



i dx 

-{di + d2)-E)g'^^\xi;x2)+t6{xi)e'^^\x2) = 



1 

--^ + f/5(x) + e^-E)e«(x)+t^7«(0;x) = 0, 
1 ax / 

(^kd, + d2)-E^g^'\x,,X2) = 0. (5) 

We construct the eigenfunctions g^^\xi, X2), g^^\xi;x2) and g^'^\xi,X2) by imposing the 
conditions that, in the region xi,X2 < 0, they are free-electronic plane waves. The eigen- 
function g^^\xi,X2) is discontinuous at xi = and X2 = 0, g^^\xi;x2) at xi = 0, and 
g(o,i)^^^ at X = 0. The value of the functions at the discontinuous point cannot be deter- 
mined by Eqs. (El). We then set g^°^{x,0) = {g^^\x,0+) + g^^\x,0-))/2 and so on. The 
function g^'^\xi,X2) should be a free-electron eigenfunction. The eigenfunctions with the 
energy eigenvalue E = ki + ^2, (fci, k2 € M) are then given as follows: 

2^(°) (xi, X2) = ^ sgn(Q) {gk, {xgjgk, (xq^ ) + uZi2(xq,q Je'-^^«2 ^(xq J) , 
Q 

e(°) (x) =gk, (x)efe2 -gk2 {x)ek^ + ^Zi2{-x)e'^'', 
g'^^\xi; X2) = gki{xi)hk2{x2) -uXi{xi2)e'^''^9{xi) , 

e«(x)=efc,/ifc,(x) + ^Xi(-x)e'^^ 
it 

2^(^)(xi,X2)= ^sgn(g)/ifcj(xQj/ifc2(xQ2), (6) 
Q 

where Q = (Qi, Q2) is a permutation of (1, 2), Xij = Xi — xj, u = 2U /(2 + iU) and 

= ih - k,)ek,ek^e'^''-'''y-'9{-x), 

X,(x) = ^e,.e^(--^*"^)^^(-x). (7) 
V 27r 



The wave-number set {ki,k2} in each of the eigenfunctions 5'^°''(xi, X2) and 5'^^-'(xi; X2) is 
not conserved during the scattering; the plane wave with {ki,k2} is partially scattered to 
that with {ed — iP, E — + iP} in the region xi, X2 > 0. In this sense, they are not the 



Bethe eigenfunctions 13l . llSl . |23| |. We have found similar eigenfunctions in the Anderson 
model [24 1 . 



The second term of each of the first four eigenfunctions ([6]) comes from the Coulomb 
interaction. The imaginary part of the wave numbers, iP, indicates the appearance of a 
two-body bound state e~^^^^^^. The interaction is a necessary condition of the appearance 
of the bound state and the strength of binding is determined by the transfer integral t. A 
similar two-photon bound state has been found in a one-dimensional waveguide coupled to 
a two- level system I25I, where the bound state has been obtained through an "S-matrix" 
acting on the Hilbert space of free two photons and the eigenstate including the bound state 
has not been constructed. 

We obtain two-electron eigenstates by inserting the eigenfunctions (E]) into the form (I3l); 
we denote them by |2, n; /^i, ^2), {n = 0,1,2). We notice that, by exchanging ki and /c2 
in \2,l;ki,k2), we have another eigenstate \2,1] k2, ki) with the same energy. The four 
eigenstates satisfy the orthonormal relations in the limit L 00: 

{2,n;ki,k2\2,n;k[,k'^) = S{ki-k[)5{k2-k'^)-S{ki-k'^)S{k2-k[), (n = 0,2) 

(2, 1; fci, A;2|2, 1; k[, k'^) = 5{ki-k[)6{k2-k'^). (8) 

In principle, we can construct eigenstates for a few electrons. For example, the three- 
electron eigenfunctions in the sector with No = are given by 

?>\g^'^\xi,X2,X'i) = ^sgn{P)gkp^{xi)gkp^{x2)gkp^ix3) 
p 

+ |5^sgn(PQ)^7.,^(xQjZp,P3(a;Q,Q3)e'('=-.+'=-3)-Q3^(a;Qj 
P,Q 

-^^sgn(Pg)V^(a;QjZp,P3(xQ,Q3)e^('=-^+'=-3)-Q3^(a;Q3Qj^(xQ,Qj^(xQj, 

2\e^°\xi,X2) = ^sgn{P)gkp^{xi)gkp^{x2)ekp^ 
p 

P,R 
P,R 

-^$^sgn(Pi?)/^fc,^(xpjZp,P3(-a;pJe^('=-^+'=-3)-H.^(a;^^^J^(a;^J. (9) 

P,R 

Here P = (Pi, P2, P3) and Q = (Qi, Q2, Qs) are permutations of (1, 2, 3) and R = (Pi, P2) is 
that of (1, 2). The third term of the eigenfunction g^^\xi, X2, 0:3) indicates a new three-body 
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bound state. The eigenstates in other sectors with No = 1,2,3 are constructed in similar 
ways. 

Now we construct a scattering eigenstate by taking a hnear combination of the four 
two-electron eigenstates as 

\ki, k2) = A\2, 0; h, k2) + Si|2, 1; fci, fca) - S2I2, 1; fca, k^) + C|2, 2; k^, fcs). (10) 

Going from the eigenfunctions in terms of the even and odd parts back to the ones in terms 
of the leads 1 and 2, we have f^^^'^\xi,X2) = (0|ci/2(x2)ci/2(a;i)|A;i, /C2) and f^^\xi]X2) = 
(0|c2(x2)ci(xi)|A;i, ^2). By choosing A = Bi = B2 = C = 1/2 in Eq. (ITU]) , we obtain the 
scattering state which contains an incoming two-electron plane wave only in the lead 1, i.e., 
/•^^^(xi; X2) = /^"^^xi, X2) = for xi, X2 < 0, which is depicted in Fig. [H In the same way, 
by choosing A = —Bi = B2 = —C = 1/2, we obtain the scattering state which contains 
an incoming one-electron plane wave in each lead, i.e., /'•'^'^-'(xi, 0:2) = for Xi,X2 < 0. We 
denote the former/latter scattering state by \ki,k2)±. Each scattering state is shown to be 
a solution of the LS equation whose incident state is a free-electron plane-wave state, where 
the incident state means an eigenstate of the Hamiltonian ([T]) with t = 0. On the other 
hand, the scattering state constructed from the Bethe eigenstates [l8|, l23| is interpreted as 
the solution associated with an incident state that depends on the parameter U. We remark 
that the scattering states are also constructed from a superposition of an infinite number of 



the degenerate Bethe eigenstates [23 1. 



We use the two-electron scattering states to calculate the quantum-mechanical expecta- 
tion value of the current / in Eq. ([2]). The expectation value with respect to the scattering 
state \ki,k2)±, {ki < ^2) is calculated as 

{ki,k2\I\ki,k2)± 27r. ^,rf, \\ , t a 1 \ 

{k,M\Kk2)^ =-(/0(fc0±/0(fc.))+-^/^(fe„fc2), 

Io{k) = ^lm(efc). 



I±{k,h) = ^ {Re{ef,)lm{uel)±Re{ek)lm{uel)) , (11) 



where L = 27r(5(0) is the length of the system. The first term of order gives the current 
of non-interacting electrons. The correction term of order containing I±{ki,k2) is due 
to the two-body bound state. 

We find that, in the limit L,N 00, the correction term in Eqs. (fTTj) contributes to 
the current. In the spirit of the Landauer formula, we assume that electrons are completely 



thermalized in each reservoir before returning to the system. We speculate from the result 
of = 2 that, for general A^, similar n-body bound states, {1 < n ^ N) contribute to 
the term of order L~" in the expectation value. We assume that the contribution from 
the two-body bound state is given by the function I±{k, h) in Eqs. flTTl) . Let \k) be an A^- 
electron scattering state with an incoming A^^f^-electron plane wave characterized by distinct 
wave-numbers {fcf } in the lead a. The speculated form of the expectation value is 

' ' ' i=l i=l 

i<j i,j i<j 

We have verified this for = 3. We neglect the terms of order higher than in the 
expansion 16|. By taking the limit L,Na oo, the sum {2tt/L) should be replaced 

by the integral on k with the zero-temperature Fermi distribution fa{k) = 9{iia~k). For 
/ii/2 = ±^/2, the average current is then given by 

l-V/2 i.V/2 .V/2 i-V/2 i-V/2 

(/) = / dkIo{k) + -[ dk dh ~ dk dh)l+{k,h) 
J-v/2 2vJ_^ J_jy J_p^ J 

CV/2 l-V/2 

+ dk dhl_{k,h), (12) 



J -A J -A 

where —A is the low-energy cut-off. We have 



J = 2(A+j+)j2 + (j--Ji)log- 



J' = 2(A+j+)ji+(j2 + ^log^^) log- 



where e± = (e^ ± V/2)/t^, ca = (e^ + A)/P, A = 2(A/? - arctan(eA)), j± = arctan(e+) ± 
arctan(e_) and js = e^^*/(e^ + 1) ~ "/(e^ + 1), (s = 1,2). The current includes higher- 



order terms in U and, at = 0, agrees with the perturbative result ^0|] in the first order in 
U. The linear divergence in A — oo is due to the linearized dispersion relation in Eq. ([T]). 
In Fig. [2], we plot the current- voltage characteristics at = by setting A = V^. The regime 



of negative differential conductance appears for large U |l9l. l21|. 

In summary, through the Landauer formula, we have studied the nonequilibrium current 
in an open quantum dot system by using exact scattering eigenstates. We have found that 
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FIG. 2: A current-voltage characteristics of the average current for = and [/ = 0, 2, 4, 6, 8, 
10. 

the effect of the interaction appears through the many-body bound states. By taking the two- 
body bound state into account, we have calculated the average current, which agrees with 
the perturbative result [20i| at = and has a behavior similar to the other results 19, 2l| . 



In order to compare our result, including the case 7^ 0, with the result in Ref. 20| precisely, 
we need to consider contributions from other many-body bound states in Eq. f|T3l) . because 



they may include first-order terms of U . They would enable us to regu 



divergences in Eq. ( fT3l) with the renormalization-group technique 19 



arize the logarithmic 



26|. 
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